Summary:
Added support for evolving binaries (no higher order multiples) with COSMIC. The objects are generated normally in SPISEA and then are evolved using COSMIC. The evolved objects are then assigned atmospheres via an interpolated atmosphere grid over a new merged atmosphere class (see below).

Details:

1. Added a new evolution object which ingests SPISEA objects into COSMIC and puts the outputs into SPISEA tables. Has the option to keep the COSMIC tables. For disrupted binaries, the user has the option to keep the disrupted companions which will be added to the primary table or to drop them. Kicks from COSMIC are also included. (Note, we use the sum of delta values of the two rows in kick_info table to accommodate both of the objects going SNe while remaining bound).
2. Added a new IsochronePhotExternalEvolution object (subclass of IsochronePhot) which creates the atmosphere grid rather than a normal isochrone. A new atmosphere grid is created with new metallicity, AKs, or distance. The grid of magnitudes that correspond with each set of parameters is interpolated with a stellar radius of 1 which is then rescaled for each individual object later in ResolvedCluster.
3. Added merged atmosphere model including black body supplement in ranges our current atmosphere models do not support (get_merged_atmosphere_w_bb_supplement).
4. Added support functions to get the ranges of all the atmosphere grids.
5. Added support in ResolvedCluster for external_evol objects which are determined by an attribute in the isochrone object. In these cases, the stars generated are passed to COSMIC to be evolved rather than using an interpolated isochrone. Note, when using external evolution, there is no IFMR used. The objects are evolved in COSMIC.
6. Tutorial here: docs/Cluster_w_COSMIC.ipynb
7. Added basic unit and integration tests with support from cursor.

Note this is already merged with the BD branch (here https://github.com/caitlinbegbie/SPISEA/tree/feature/BDs)